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Abstract 


The envelope theory, also known as the auxiliary field method, is a simple technique to compute 
approximate solutions of Hamiltonians for N identical particles in D dimensions. The quality of 
the approximate eigenvalues can be improved by adding a free parameter in the characteristic 
global quantum number of the solutions. A method is proposed to determine the value of this 
parameter by comparing the eigenvalues computed with the envelope theory to the corresponding 
ones computed with a V-body generalization of the dominantly orbital state method. The accuracy 
of the procedure is tested with several systems. 
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I. INTRODUCTION 


The envelope theory (ET), also known as the auxiliary field method, is a procedure to 
obtain approximate solutions of A^-body Hamiltonians in quantum mechanics [l-4]. In the 
most favorable cases, an approximate eigenvalue is an analytical lower or upper bound. 
The method relies on the replacement of the Hamiltonian H under study by an auxiliary 
Hamiltonian H which is solvable, the eigenvalues of H being optimized to be as close as 
possible to those of H. Recently, the accuracy of the ET for eigenvalues and eigenvectors 


las been tested by computing the ground state of various systems containing up to 8 


a. 


Dosons 
g and 


This comparison was possible thanks to accurate numerical results published in 
obtained by using explicitly correlated Gaussian basis [^. 

In several cases, it has been shown that the approximate eigenvalues computed with the 
ET can be improved by modifying the structure of the characteristic global quantum number 
of the method j^, 8]. The price to pay is generally the loss of the variational character of 


the approximation. In a modification of the global quantum number is proposed and 
tested by introducing only one supplementary parameter cp. This procedure relies on the 
result obtained in j^, where a universal effective quantum number for centrally symmetric 
2-body systems is proposed. The new parameter (j) can be determined by a fit on analytical 


B- 


results for = 2 or on accurate numerical results 

In this paper, we propose a method to compute the value of 0, independently of the 
knowledge of analytical or numerical previous calculations. The idea is to compare the 
eigenvalues computed with the ET to the corresponding ones computed with a A^-body 
generalization of the dominantly orbital state (DOS) method. In this last method, an ap¬ 
proximate solutions is found by quantizing the radial motion around a semiclassical solution 
for a circular motion. Developed at the origin for 2-body systems , it has been extended 
to 3-body systems jn |. 

The paper is organized as follows. The ET method is briefly recalled in Sec. [TTl The 
extension of the DOS method for A^-body systems is presented in Sec. IIIII In Sec. HVl the 
computation of the parameter 0 is described. The relevance of this parameter is tested with 
various systems in Sec. |Vl Since the behavior of observables is not really predictable when 


0 is varied 
given. 


B. 


we do not consider eigenvectors here. Some concluding remarks are finally 
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II. THE ENVELOPE THEORY 


We consider the following general Hamiltonian, in a D dimensional space {D >2), for N 
identical particles 

N N N 


H = J2T{\Pi\) + J2^{\ri-R\)+ Y1 V(\ri-r,\). 

i<j=l 


( 1 ) 


i=l 


i=l 


T is a kinetic energy, U a one-body interaction and V a two-body potential {h = c = 1). The 
center of mass motion is removed Pi = 0) ^ is the center of mass 

position. It is shown in j4| that, in the framework of the ET, an approximate eigenvalne E 
is given by the following set of eqnations for a completely (anti)symmetrized state: 


VCn J ’ 


E = NT{po) + iV f/ (^) + H 

roPo = Q, 

Np„r(p„) = r„f/' (^) + Vc^r„V' (-^) . 


\VCn J 

where W\x) = dW/dx, Cn = N{N — l)/2 is the nnmber of particle pairs, and 

N-l ^ 


( 2 ) 

( 3 ) 

( 4 ) 


Q — -f- li) + {N — 1)-— 


( 5 ) 


2 = 1 


is a global qnantnm nnmber. This corresponds to iV — 1 identical harmonic oscillators, 
implying a strong degeneracy. Following the forms of T, U and V, the approximate valne 
E can have a variational character jd]. 

The method is easy to implement since it rednces to hnd the solntion of a transcendental 
eqnation. Eqnations dSM]) have a nice semiclassical interpretation jd], bnt the ET yields trne 
qnantnm resnlts: approximate eigenvalnes and eigenvectors. In this paper, we focns on the 
eigenvalnes. Information abont the corresponding eigenvectors can be fonnd in j^, . Let 

us simply mention the following results giving a clear interpretation of the quantities Pq and 

ro Isl: 


N 


K=n{'Ep‘)’ 


. i=l 
N 


N 


^0 = 




( 6 ) 

( 7 ) 


. 2=1 


\i<j = l 
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III. THE iV-BODY DOMINANTLY ORBITAL STATE METHOD 


In the DOS method for 2-body systems, an approximate solution is found by quantizing 
the radial motion around a semiclassical solution for a circular motion [^. In 111, this 
method has been extended to 3-body systems, but it is mentioned that a generalization to 
A-body seems difficult. Nevertheless, this appears hnally possible starting from ([2111]), since 
a semiclassical interpretation of the ET is established in terms of circular motions of N 
particles in a fully symmetrical conhguration j4|. 

Let us consider a pure circular motion for the N particles with the energy Eq. We can 
rewrite (|3]) as roPo = -^ 4], where the total angular momentum A will be specihed later. 
But, it is expected that A S> 1 in order to justify the semiclassical approximation. Eq can 
be computed from the system (12}|1|) with Q replaced by A. Now, let us determine a solution 

ion around this circular 
111, we build from ([2]) a 


with energy E = Eq + AE for a small collective radial perturba 
solution. Following a similar procedure as the one developed in 
Hamiltonian AH ruling this radial motion by the replacements, 

ro ^ rQ + Ar, 


Po = — ^ \IPr + 
ro 


A2 


( 8 ) 

(ro + Ar)2’ 

where Pr and Ar are conjugate variables. Assuming that (Ar) tq and (p^) "C po, a power 
expansion of AH, limited to the lowest non-vanishing orders, gives 

( 10 ) 


AH ^ —pi + —Ar^, 

2/" 2 ’ 


with 


A 


/i = 

k = 


AroT'(A/ro)’ 

^ (2roT'(A/ro) + AT"(A/ro)) + ^M^vq/N) + H"(ro/^). 

Tn A 


( 11 ) 

( 12 ) 


The terms in Ar are cancelled by (01). The perturbed energy AE is then simply given by 


AE = Ah' with A = 


and u = 


1 3 
2 ’ 2 ’ 


(13) 


It is expected that i/ -C A in order that AE Eq. Explicitly, 


A2 = 


2N^ 




T'2 ( ^ 


+ 


N^X 


A 




r —] T" — 


A 


ro 


A 


A 


ro 


^T' ( — ) U" i^] + 


N 


Nro 

A 


T' ( — 




^^0 


V Cjy 


(14) 
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with ro being a solution of the system (EJHj) where Q is replaced by A. To produce relevant 
results, quantum numbers A and u must be specihed. This work is done for = 3 in Q. 
The general case is discussed in the following section. Let us recall that the results given 
by this method are expected to be relevant for completely (anti)symmetrized states [^|, 
provided i/ -C A and A 3> 1. 


IV. IMPROVEMENT OF THE GLOBAL QUANTUM NUMBER 

In { 9 ], it is shown that an effective quantum number q that determines with high accuracy 
the level ordering of centrally symmetric 2-body systems has the following form (for practical 


purposes, our dehnition of 0 differs from the one used in 


m 


, ^-2 

q = (j)(n+- \ +/ + _—, 


(15) 


where the number 0 depends on the system. This structure comes from the separation of 
angular and radial motions. The DOS method relies on a similar separation for N = 2 and 
3 [ 10 , lul • For V-body systems, we can try 


N-l 


jj = + (A^- 1)-, 


2=1 

7V-1 

X = Y,k + {N-l) 
2=1 


2 

D-2 


(16) 

( 17 ) 


This is the case in [m| with V = 3, where the quantum numbers z//\/3 and A/\/3 are used 
due to a particular scaling of the conjugate variables. 

It is shown in j^, that the energies obtained by the ET can be signihcantly improved by 
replacing the quantum number Q in the system (EM]) by the new quantum number 


= (j)u + X, 


(18) 


which partly breaks the strong degeneracy of Q. This is an obvious generalization of 
but for V — 1 degrees of freedom. Such an extension from 2 to V particles is expected to be 
relevant because of the existence, in the framework of the ET, of relations linking the energy 
of a V-body system to the energy of the corresponding 2-body system by a rescaling of the 


global quantum numbers 


12| . The genuine ET, with its possible variational solutions, is 
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then recovered with 0 = 2. For other values of 0, the variational character of the solution 
cannot be guaranteed. 

If the parameter 0 is assumed constant, it can be determined by a £t on a single known 
accurate solution. This procedure is used in for some specihc systems. It should be 
preferable to have procedure allowing the computation of 0 without the prior knowledge on 
some solutions of the system. This is possible by combining the information coming from 
the ET and the DOS method. 

Let us write = A(1 + e) with e = 0z//A. Ifi/ = e = 0, the solution Eq of the system ([2]- 
0 ]) is the unperturbed solution of the DOS method. We can now search for a solution when 
e -C 1, by setting 


ro ^ ro + Ar, 

A A(l + e) 

Po — - 1 — . . 

tq tq + Ar 


(19) 

( 20 ) 


An expansion to the non-vanishing lowest orders of (jl]) yields a relation between e and Ar. 
Using this link in an expansion to the non-vanishing lowest orders of ([2]) hnally gives 


AE = Be, 


( 21 ) 


with B = Bn/Bd and 

NX 


B„ = 




2T' ( —^ + —T''(— 

roj ro \ro 




ro 


+xr (- 

To 


U"{^]+NV" 


ro 




r, 


0 


ro 


r, 


0 


ro 




N 


VCn 

Equations flT^ and correspond to the same expansion z/ A. So we can write 

, AA 


B 


( 22 ) 


(23) 


(24) 


The complete expression for 0 is quite complicated. But, the computation of 0 is straight¬ 
forward once ro is determined. In some cases, the hnal expression can be very simple, as 
shown with the examples below. 

The procedure to compute a numerical solution is the following: 
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1. Fix the global quantum numbers v and A; 


2 . Compute the solution vq of the system ([3}11]) with Q replaced by A; 

3. Compute (p with fl21D : 

4. Compute = 0 1 / + A; 

5. Compute the solution E of the system ([2}11]) with Q replaced by Q(j). 

If ro(Q) can be determined analytically, 'ro(A) is used for step 2 and is used for step 

5. 

With the replacement of Q by Qfp in the ET equations, we can hope that some relevant 
information about the relative influence of radial and orbital motions are now captured in 
the ET. Nevertheless, (j) is determined in a particular limit: z/ A and A S> 1. It is not 
clear that the dehnition of can yield good results for other ranges of quantum numbers. 
This will be tested in the next section. 

V. RESULTS 

In this section, we consider various systems of N identical particles with the same mass 
m, for which analytical ET solutions can be found. The ET formalism developed is valid for 
bosons or fermions with any values of D, but comparisons with numerical results will only 
be performed for boson-like particles in the D = 3 space, as in [^. 

A. Power-law 2-body potential for nonrelativistic particles 

In the hrst example, nonrelativistic particles interact via a 2-body power-law potential 



(25) 


with a > 0 and 0 b > —2. The ET gives the following upper (lower) bound for the energy 
E iib < 2 {b> 2) 



( 26 ) 
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The computation of 0 gives a remarkably simple result 


0 = Vb + 2 . 


( 27 ) 


When b = 2, the exact result 0 = 2 is recovered. 

Let us dehne the ratio R = lim,^>A E/ hmA>j. E. Formulas (l26|) and (l27|) give 


i?ET — Ci{b) j 


l/\ 2b/{b+2) 


with Ci{b) = {b + 2f/^^+‘^\ 


(28) 


The Hamiltonian (12^ has been studied for N = 2, D = 3 and 6 > 0 by a Bohr-Sommerfeld 

n 

quantization (BSQ) approach in |13|. From formulas (12) and (14) of this reference, one can 
get with our notation 



( 29 ) 


where B{x,y) denotes the beta function. The dependence on quantum numbers of these 


two ratios is the same since v = ni and \ = li asymptotically when N = 2. The coefficients 


Clip) and C 2 {b) are different. But, for b G [0,2.5], the relative difference 5{b) = |C'i(6) — 
C 2 {b)\/{Ci{b) + C 2 {b)) does not exceed 1.6%, with 5(0) = 5(2) = 0. For larger values of b, the 
situation is less favorable, since limf,^oo (Fi(6) = oo while hm6_j.oo ^2(5) = vr^. Nevertheless, 
the agreement between the two approaches can be considered as quite good, since 0 is 
computed with the assumption that X ^ u, not u :$> X. 

We have checked that 0 = \/2 if V{r) = a ln(r/&). This corresponds to 6 = 0 in (l27|l . as 
expected for a logarithmic potential. 

B. Power-law 1-body interaction for ultrarelativistic particles 

In the second case, ultrarelativistic particles (m = 0) interact via a 1-body power-law 
potential 


r(p) = p, t/(s) = as*, V'(r) = 0, 


( 30 ) 


with a > 0 and 5 > 0. The FT gives the following upper bound for the mass F if 5 < 2 (if 
b>2, E has no determined variational character) 


ER-^(NahQ’i) 




(31) 
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Again, a very simple result is obtained 


0 = V6 + 1. (32) 

We have checked that (j) = 1 ii U{s) = a ln(s/6). This corresponds to 6 = 0 in fl3^ . as 
expected for a logarithmic potential. 


C. Weakly interacting particles 


Let us now consider nonrelativistic particles interacting via a soft 2-body Gaussian po¬ 
tential 


r(p) = 1^, U(s) = 0, V{r) = 


With the dehnition 


= -VITUivV 


{N - 1)R./2^q 

the resolution of the ET equations gives the following upper bound for the energy E 


E = 


A'(iV - 1) „ ,, l + 2 tr„(F(Q^)) 


VoYiQtY 


(33) 


(34) 


(36) 


2 W„(Y(Qt)f ' 

where the multivalued Lambert function W{z) is the inverse function of z Q- Wo(z) is 
the branch dehned for z > — 1/e. The computation of (j) gives 


.#. = 2^1 + M-'o(r(A)). 


(36) 


When R — )■ oo, the interaction reduces to a constant plus a harmonic potential, V{r) k, 
—Vo(l ~ IR?\ in the physical range of the eigenstates. It can be checked that E —)• 
—Gat Vq + ■\/2 iV Vo/ m R^ and 0 —)• 2, as expected. 


The good quality of the ET results for this kind of systems has been checked in [5 


for the 


ground state of bosons, thanks to the accurate numerical results obtained in 0,ll^. With 
the parameters used in [^, fl36|l gives values of 0 varying from 1.85 to 1.95, when N runs 
from 2 to 20. This is in good agreement with the constant value 0 = 1.82 htted in j^. 


D. Self-gravitating particles 


The following Hamiltonian, 




( 37 ) 


9 










is a particular case of (125|) with h = —1 and 0 = 1. For a self-gravitating system, g = im?G 
where G is the gravitational constant. Thanks to the accurate numerical results obtained in 
it has been shown in that the ET results for the ground state of systems composed 
of up to 8 bosons can be largely improved by a value of 0 around 1. The value 0=1 also 
allows a good agreement between the ET prediction and a lower bound for the ground state 
of N bosons in the D = 3 space H- 


E. Confined particles 


Now, we consider particles confined by a harmonic oscillator potential with a pairwise 
repulsive Coulomb interaction 


Up) = G’ = -■ 

2m 2 r 


With the definition 


F(Z) = 


216/3 


u 


3 m/^{N-ly 
the ET method gives the following lower bound 

, , ,, 1 /s 


2/3 




E = 


25/3 


(G4y{Qi,)f + G-(y («*))"') ■ 


(38) 


(39) 


(40) 


where G±{Y) is the only positive root of the quartic equation 4a;^±8a; —3E = 0 with E > 0 
h For the bosonic ground state, the value 3uj/2 must be added to fl40|) if the confinement 
is ruled by ’"i instead of — Rf. The computation of 0 gives 


^24/M^ + l. 


E(A) 


(41) 


One can check that E = 00 with 0 = 2 when g = 0, which is the exact solution in this 
case. 

The comparison with the accurate numerical results obtained in jo] for the ground state 
of systems composed of up to 8 bosons with the parameter m = 1, a; = 0.5 and g = I 
shows that the ET gives good lower bounds. Better results for the ET can be obtained 
by increasing, with the number of particles, the value of 0 above 2. This is achieved by 
formula flTT|) . but too much to give an improvement of the results. This is illustrated on 
Fig.[Il 

^ G±(y) =T|\A^+|v'4(F(y))-i/2_y(y) with V{Y) = (2 + V4 + W)- F (2 + V4 + 
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Figure 1. Log-log plot of energies, E, of N confined bosons: exact results (circle) 
for (p = 2 (diamond), ET results for </> given by (j4ip (dashed line). 


li, 


ET results 


F. Large-A^ baryons 


In a particular large-iV limit of QCD {N is the number of colors), light baryons are com¬ 


posed of N ultrarelativistic quarks u or d {m = 0 


In this case, the corresponding Hamiltonian is 17| 


in the fundamental color representation. 


T(p) = p, U{s) = ks, V (r) = — 


9 


(42) 


Quarks are fermions but they can be treated as bosons since the color part of the wave- 
function is completely antisymmetrical. With the ET, the following upper bound is obtained 
for the mass E of the baryons 

' 3/2 


E = \/4^ 

The computation of 0 gives 




9- 


= \ 2 - 


^N{N-irg 

V2X 


(43) 


(44) 


These results with g = 0 are in agreement with those of Sec. IV Bl with 6 = 1. 

In order to test the relevance of for the ground state and some excited states, we 
consider the 3-body system studied in [s]. For the genuine ET (0 = 2), the mean relative 
error A computed on the 16 eigenmasses listed in Table [H is 15.1%. With the value 0 
predicted by (011), A reduces to 4.7%. More precisely, the mean error on the 4 states of 
Tabled] with /i -|- /2 = 0 is 9.1%, while it is only 3.2% on the 12 remaining states. This could 
be interpreted as a trace of the condition of validity of the DOS method (A 1). 
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Lets us remark that better global agreements can be obtained with constant values of 0. 
With 0 = 1.35, the ET yields the exact value for the ground states and A drops to 3.1%. 
Finally, with 0 = 1.23, the minimal value of 2.4% is reached for A. Obviously, it is necessary 
to known some exact solutions of the system to compute these values of 0. From Table [U it 
is also clear that the variational character of the ET solutions is lost when 0 % 2. 


VI. CONCLUSION 


The envelope theory is a powerful method to compute eigenvalues for quite general N- 
body systems with identical particles in D dimensions The method is easy to implement 
and can produce analytical upper or lower bounds in favorable situations, as those presented 
in this paper. Its purpose is not to compete with accurate methods such as the explicitly 
correlated Gaussian basis or the harmonic oscillator basis [^, but to produce without 
great pain reliable estimations for the energy. 

It has been shown in that the eigenvalues computed with the envelope theory can be 
quite close to the exact results, and that they can be improved by the introduction of a 
single parameter 0 in the characteristic global quantum number of the method. This pro¬ 
cedure has been inspired from Qj, where a universal effective quantum number for centrally 
symmetric 2-body systems is proposed. The price to pay is that the variational character 
of the approximation cannot then be guaranteed. This parameter 0 can be determined, for 
instance, by the knowledge of an exact solution for iV = 2 or by a ht on numerical results. 

In this paper, a method to compute the value of 0, independently of the knowledge of 
analytical or numerical previous calculations, is proposed by using, in the envelope theory, 
information coming from a A-body generalization of the dominantly orbital state method 
I^. 11]. A procedure to compute 0 is given in the general case. Though, the complete 
formula is quite complicated, it can reduce to very simple forms, as for the cases presented 
in this paper. 

From the way 0 is computed, it is expected that the improvement of the eigenvalues is 
good only when the global angular momentum is large and the global radial excitation is 
small. But the situation is not so clear, since even ground states can be improved. With a 
value 0 % 2, the variational character of the approximation cannot be guaranteed. It is then 
not possible to locate with certainty the approximation with respect to the exact energy. 
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Nevertheless, the artihcial strong degeneracy inherent to the envelope theory is broken, and 
relevant insights about the relations between orbital and radial excitations can be obtained. 

Since iV-body problems are always difficult and heavy to solve accurately, the envelope 
theory, supplemented by the computation of 0, can be used as a guide for the study of these 
complicated systems. It could be interesting to test this method with other systems in order 
to better specify the domain of validity of this procedure. Accurate numerical results for 
various systems, with different values of N and D, are welcome. 
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Table I. Some eigenmasses in GeV of the Hamiltonian (14211 for D = 3, N = 3, k = 0.2 GeV and 

n 

g = with as = 0.4 [^. Accurate results obtained from an expansion in a harmonic oscillator 


18l | are given in the column “Exact”. Results computed with the ET for four values of cj) are 


listed in the four last columns, with the associated mean relative errors A. The sums ni + n 2 and 
li + I 2 are the quantum numbers of the main component of the corresponding genuine eigenstate 
in the harmonic oscillator basis. These numbers are used for the computation of 


ET 


ni + 

n2 h 

+12 Exact 

(/) = 2 

dm) 

(j) = 1.35 (j) = 1.23 

0 

0 

2.128 

2.468 

1.945 

2.128 

2.060 

0 

1 

2.606 

2.914 

2.582 

2.633 

2.578 

1 

0 

2.739 

3.300 

2.504 

2.788 

2.682 

0 

2 

2.959 

3.300 

3.035 

3.055 

3.007 

1 

1 

3.125 

3.646 

3.106 

3.189 

3.098 

0 

3 

3.299 

3.646 

3.418 

3.425 

3.383 

2 

0 

3.260 

3.961 

2.960 

3.318 

3.186 

1 

2 

3.422 

3.961 

3.512 

3.546 

3.463 

0 

4 

3.581 

3.961 

3.758 

3.759 

3.721 

2 

1 

3.584 

4.253 

3.553 

3.662 

3.542 

1 

3 

3.716 

4.253 

3.857 

3.869 

3.794 

0 

5 

3.861 

4.253 

4.068 

4.066 

4.030 

3 

0 

3.721 

4.527 

3.354 

3.775 

3.619 

2 

2 

3.838 

4.527 

3.932 

3.976 

3.866 

1 

4 

3.966 

4.527 

4.166 

4.168 

4.098 

0 

6 

4.103 

4.527 

4.356 

4.351 

4.318 

A 



15.1% 

4.7% 

3.1% 

2.4% 
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